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Abstract 

Using classical electrodynamics we determine the angular dependence of the 
light intensities radiated in second and third harmonic generation by spherical 
metal clusters. Forward and backward scattering is analyzed in detail. Also 
resonance effects in the integrated intensities are studied. Our work treats 
the case of intermediate cluster sizes. Thus it completes the scattering theory 
fo spherical clusters between Rayleigh-type analysis for small spheres and 
geometrical optics for spheres much larger than the wavelength for nonlinear 
optics. Since the particle size sensitivity of Mie-scattering is increased by 
nonlinearity, the results can be used to extract sizes of small particles from 
nonlinear optics. 
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I. INTRODUCTION 



In this paper we investigate the nonhnear interaction of hght with small metal particles. 
In nonlinear optics one expects a more pronounced sensitivity of the radiated yield to the 
light-wavelength and cluster-size, due to effects by nonlinear sources. 

The linear interaction of light with spherical objects shows certain changes depending on 
the size of the sphere. The limiting cases are Rayleigh- scattering for spheres small compared 
with the wavelength and reflection by spheres much larger than the wavelength (the limit 
of geometrical optics). If the size of the sphere and the wavelength of light have comparable 
magnitudes the intensities scattered by spherical metal clusters are strongly enhanced. Mie 

first calculated this classical effect and found that most of the incident light intensity is 
scattered in the forward direction, because of interference of waves originating from the front 
and the back of the sphere. This so called Mie-effect is increased with increasing sphere size. 
Bohren and Huffmann pointed out that the extinction efficiency Qext defined as the sum 
of scattering and absorption efficiencies {Qext = Qsc + Qabs) as a function of the sphere-radius 
are characterized by the interference- and ripple structure. The interference structure gives 
rise to regular interferences of incident and forward scattered light dominating the envelope 
of Qext^ see Fig. (|1|). Superimposed on this is sharp and highly irregular ripple structure 
originating from resonant electromagnetic surface modes of the sphere. In general, the 
resonances depend on the size parameter ka comparing the radius of the sphere a and the 
wavelength expressed by the absolute value of the wave vector k of the incident light. Thus 
resonances also occur by varying the wavelength of the incident light. For small spheres only 
one peak or resonance appears at ujp/ with ujp the plasma frequency and the 1 / v^-factor 
caused by the shape of the sphere. This is the dipole-term. The larger the size of the sphere 
the more resonances around the cUpZ-y/S-peak appear (if ka > 1) reflecting the surface-mode 
resonances or higher multipoles, see 0. 

The nonlinear interaction of light with small particles was recently studied by Ostling 
et al. 1^ using a simplified model to obtain the intensities in second and third harmonic 



2 



generation, SHG and THG, respectively. The calculations yield an enhancement of the 
total intensities by a factor 5000 in SHG and 200000 in THG in comparison to a plane metal 
surface for spheres with sizes a > Xp/2 in a wide frequency range. Here, Xp is the wavelength 
corresponding to ujp using the relation c = \ ■ Lojl-n with c the speed of light in vacuum. 
This holds especially if the diameter of the sphere roughly equals \p. 

The correlation between the scattering behavior and the properties of the sphere enables 
us to obtain the size and complex refractive index of the sphere and in the case of an 
ensemble of spheres even the size distribution from the scattered intensities. It is a goal 
of our investigation to show that this is improved by higher harmonics. In this paper 
we extend the theory of ref. ^ to obtain the angular dependence of the radiated second 
and third harmonic intensities, compare the linear and nonlinear results examining their 
characteristic features in higher harmonics. In particular the angular resolved intensities, 
the degree of polarization of the scattered SHG and THG yield are investigated. In addition, 
we thoroughly study parameter dependences of the ratio of forward and backward scattering 
and analyze the origin of the resonance structures. 

Our main results are an enhanced size sensitivity of the higher harmonics compared to 
linear Mze-scattering, thus pronouncing the mentioned characteristic features of the linear 
Mie-effect in higher harmonic scattering. Similarities between the linear scattering and THG 
reflect that the scattering is dominated by the same multipoles in both cases. In contrast, 
different multipoles contribute to SHG resulting in distinct changes of the angular depen- 
dence. Because of the more pronounced forward scattering enhanced backward scattering 
vanishes in the higher harmonics. 

The following section II contains the appropriate theoretical calculations. In section HI 
we present the numerical results starting with polar plots of the intensities in SHG and 
THG, plots for the degree of polarization, the ratio of forward and backward scattering and 
the integrated intensities as a function of the size of the spheres. In section IV we will give 
a summary of the results. 
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II. THEORY 



First, we use Mie theory to calculate the linear scattering and the radial component of the 
polarization a {6, 0) at the surface of the sphere. Then, we determine the radiated intensities 
in second and third harmonic by matching the electromagnetic fields and the n-th power 
of cr {6, (j)) with regard to the appropriate boundary conditions, as described by Ostling et 
al. In section II. D the characteristic quantities of the radiation are determined. These 
quantities can be compared with experiments. 



A. Linear scattering by spheres 

According to the spherical symmetry of the problem we expand the fields in the form 

(x) = ^ C (/) \a\, (l, m) // (hr) X,,^ {6, 0) 

(1) 

Therein, X;_m is a vector spherical harmonic as introduced by Jackson |^ with C (/) = 



Lm 



Tfl 1 

+ (/, m) — -y V X /; (fcir) Xi,„ {9, , 



'An {21 + l),k = uj/c and ki = ye [uj)k . The multipole coefficients (/, m) and (Z, m) 
refer to the magnetic (transverse electric) and electric (transverse magnetic) multipoles. The 
index i specifies the incident {i = inc), the scattered {i = sc) or the internal {i = in) fields. 
For incident waves of positive and negative helicity we have a^^f{l,±l) = a^'^(/,±l) = 1. 
In this paper we use a superposition of both to give linear polarization. The spherical 
Hankel functions ff^ {kr) = /i^^ (kr) and Bessel functions J*"'^'™ {kr) = ji {kr) describe the 
radial part of the field outside and inside the sphere. The magnetic field is given by the 
Maxwell-equation for harmonic fields 

B = -iuj/c-V xE. (2) 

Using the boundary conditions at the surface of the sphere 
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n X (e^^ + E*"^) = n X E^ 



(3) 



and 



n X (b^'^ + B^"^) = n X B^", 
we obtain the expansion coefficients of tlie scattered wave 
af (Z,±l) - 



(4) 



ji {kr) £ [rji {kir)] - e {uj)ji {kir) £ [rji {kr)] 



e {ui) ji {kir) £ [r/ip^ {kr)] - /i^^ {kr) £ [rji {kir)] 



a-(/,±l) = 



Ji (kr) §. [rji {kir) 


- ji {kir) £ [rji {kr)] 


ji {kir) i 


rhf^ {kr) 


- h^^ {kr) I [rji {kir)] 



(5) 



at tlie surface of tlie spliere witli radius a. From tlie continuity of tlie electrical displacement 
at the surface of a perfect conductor 



n • (D"^ + W""") = n • D^'', 



(6) 



the surface charge results as 



-iijjt 



(7) 



where n = r/ |r| and Re denots the real part. Furthermore we expand a {9,(f)) in spherical 
harmonics 



<^(^,0) = 2 S ai%Yi,^{0,<f>)e-^' + c.c. . 

l,m=±l 

The expansion coefficients result from the orthogonality of the spherical harmonics as 
1 



(8) 



47r 



C{l)iJl{l + l) , , 
\^ {ji {ka) + a- {I, ±1) /i!^) {ka)) . 



(9) 



B. Sources of the higher harmonic radiation 

In analogy to the linear case we expand the n-th power of the surface charge a in terms 
of spherical harmonics: 
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{0, <P) = lJ:<^^S^^l,m {0, <P) e-"^-* + C.C. . (10) 
l,m 

Neglecting time-independent terms we obtain the coefficients in the case of second harmonic 
generation as 



-j^ oo oo „ 
-j^ oo oo „ 

~ 2 ^ ^ ^h-l'^h-l J ^l*-2^li,-l^h,-ld^ > 



a 



1,-2 



2 

^ il=l l2 = l 



1 oo oo „ 

?o = 2 E E j Y*,Yi,,,Y,,,.,dn , (11) 



.(2) _ _ 
9 

il = lZ2=l 

and for third harmonic generation (also neglecting terms with e~^^) 



1 oo oo 

4-1 = E E «S«£-2 / Y*_,Y,,,,Yi,,.,dn , 

-j^ OO OO „ 

9 E E ^li,-l'^l2,2 / ^C3^^l.l^'2,2'^^ , 
^«i=li2=l 

-| oo oo „ 

2 E E y r,:_3y5,_,F,„_2dJi . (12) 



(3) 
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* ~ 9 

/l=l /2 = 1 

The integrals can be expressed by the 3j-symbols and yield the coupling of the multipoles. 
Because of conservation of angular momentum, only coefficients with m = 0, ±2 in SHG 
and m — ±1, ±3 in THG differ from zero. 

C. Higher harmonic radiated fields 

In the case of higher harmonics the electric and magnetic fields inside and outside the 
sphere and the sources are matched by the boundary conditions 

n . (^D°"* - D*") = 47r(7'^ {9, (f) , (13) 

and 

n X (e°"* - E^'*) = . (14) 
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As a result of spherical symmetry, only transverse magnetic waves are generated by the 
oscillating surface charge. Thus the fields in the nonlinear case are 

1 



m) 



Lm 



e (nuj) k 



(15) 



where i = out or i = in, respectively, and k = nu/c, ki = {nuj)k, //" (kr) = ji {kir) and 
jout _ ^(1) (j^^y The boundary conditions give the coefficients of the radiated field 

[rji {kir)] Txka 



4") (/,m) 



In this case we have ki = Je {nuj)k. 



d_ 

dr 



An) 



(16) 



D. Calculation of quantities characterizing the radiation 



To study the angular dependence of the scattered field we use the quantity {6, + 
\Ee {6, (f))f according to Born and Wolf where Eg {6, (p) and E^f, {6, 0) are the tangential 
components of E'^^ {6, 0) in the linear case and of E°"* {6, 0) in the nonlinear case. This defini- 
tion is equivalent to the absolute value of the radial part of the Poynting vector |n ■ (E x H) |. 
The following formulas represent l-E^ {9,(f))f and \Eg {9,(f))f in the far field approximation. 
We obtain after evaluating the m-summation in the linear case 

\E.(e.6)\^ 



\Ee{9,<P)\' 



y do sm u 

1=1 



sin 9 



d9 



■ sm 



■ cos 



for second harmonic generation 



An (21 + 1] 

1=1 



cost 



4'^ (/,0) 



\ d9 
dPl (cos 6 



2K (/) {I, 2) cos (20) 



\Eei9,<P)f 



1=2 



d9 

-dPl (cos I 



4) (/,2) 



sin^ (20) , 



(17) 



(18) 



and third harmonic generation 



\E4>{e,<p)[ 



t i^^^Af (l, 1) sin 

1=1 \ ^ 



sin 9 



+ 



dPi^ (cose) 

df 



K (l) A^E^ {1, 3) cos (30) 



(19) 



where K (l) are and m-dependent factors. 

Note that the Legendre polynomials Pf^ with m — 0,2 are identically zero for 9 = 0,7r, 
since there is no direct scattering in forward nor backward direction in the second harmonic 
case. Furthermore the 0- dependence of the hnear scattering and THG is described by the 
interval (0, tt) and by (o, | j in SHG according to the symmetries of the sine and cosine 
terms. 

We will plot the degree of polarization defined as 

III -/i 



P{9) 



(20) 



with 



Ee{9, 



+ 



(21) 



and 



\Ee{9A 



f + \E^i9,<P = 0)f 



(22) 



To measure the asymmetry of forward and backward scattering in the Mie-range we intro- 
duce the quantity 

I forw ^back 



(23) 



I forw ~l~ ^back 

which we call the "degree of Mie-asymmetry" . In the linear and Iback are the scat 



tering intensities taken at ^ = and 9 — tt respectively. As these quantities are identically 
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zero in SHG, we use for Iforw and hack the maxima of the scattering intensities along the di- 
rection of propagation of the incident wave for = 0, ^, | and 9 covering the interval (0, vr). 
In THG the angular dependence fo the radiated intensities is more complicated compared 
to the linear case and we take the maximum for = and 9 ranging from to tt. 

Finally, we calculate the angle-integrated scattered intensities. We obtain in the linear 



case 



Qic^ = ^^E7l^(l«^^(^>^)l' + Karn)|^) , (24) 

^ I'^^J i,m 



and for the n-th harmonic 



« = -^EK"^(^H' ■ (25) 



',m 



Here Q^"^ is in units of the geometric cross section of the sphere na^. In this formulation the 
optical theorem, which links the extinction efficiency to the normalized scattering amplitude 
in forward direction, has the form 

Qe.t^^—,-{\Ee{9 = 0)\'+\E^{9 = 0)f) , (26) 



with 



Qe.t = — f-2 E (2^ + 1) [af (/, m)] + Re [a^ (/, m)]} . (27) 
^{ka) i,m 
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III. NUMERICAL RESULTS 



We present numerical results for the angular dependence of the radiated intensities ob- 
tained using Eqs. (p^-(p!9|). The degree of polarization characterizing the radiation, for 
example its angular dependence. The degree of Mie- asymmetry R calculated using Eq.(p3D 
gives us the strength of the asymmetry along the direction of propagation according to the 
Mie-effect. The integrated intensities give us the resonances as a funtion of cluster size. In all 
cases we compare the linear results with the results of second and third harmonic generation 
as a function of the size parameter ka and the material-properties, referring to iron and 
nickel at a fixed optical wavelength of 617 nm. To check the numerical accuracy we compare 
the linear results with those of ref. 0. We find excellent agreement. Furthermore, we check 
the optical theorem. By determining AQ = Q^xt — 2 ■ (l-Eel^ + \E^f^ ^ ^/ n (ka)"^ as a 
function of ka we find that this is satisfied to an accuarcy better than 10^^^. 
As input parameter we take the complex refractive indices measured by Johnson and Christy 
0. In the linear case, we also use for comparison with other calculations the refractive index 
of water droplets, as given by Bohren and Huffman 0. These values are listed in Tab. I. 
The refractive index is constant in all figures unless something else is specified. Thus varying 
the size parameter ka means varying the size of the sphere. 



A. Angular dependence of the intensities 

First, we show polar plots of \E(f,f + \Egf (see Figs. H - U). In general, the shape of the 
plots is governed by the values of the coefficients in the series expansions Eqs. ([T7|)- ([l9|) . It is 
well known that in the linear series (|l^ only terms with I < ka contribute significantly. 
For / > ka the terms decrease very rapidly, whereas for I ^ ka they have comparable 
magnitudes. We restrict our calculation to / < [Max (n ■ ka, n' ■ ka) + 15] with a maximum 
value of / of 50 where n and n' are the real and imaginary parts of the complex refractive 
index. This gives satisfactory convergence of the series up to ka < 10 in all harmonics. 
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To qualify the angular dependence the range of size parameters up to 5 is sufficient. In 
this range the /-values of the dominating terms are a little smaller than ka and no terms 
with / <^ ka exist. Terms with / < ka can be very different from each other in contrast 
to terms with I -C ka which have comparable magnitudes. Thus, a pronounced transition 
range from pure dipole-scattering to pure Mie-scattering exists. Of course we cannot reach 
the transition from Mie-scattering to the optical limit of reflection, due to the numerical 
limit of a maximum value of / of 50. 

The geometry of the scattering is specified in the inset of Fig. |^ with the direction 
of propagation of the incident wave being parallel to the positive y-axis and polarization 
along the positive z-axis. Fig. ^ (a) shows Rayleigh-scatteimg according to the dipole-term 
with / = 1. The characteristic cos^ 6'-dependence appears along the x-z-plane. The other 
plots for linear optics (Figs. |^) show the well known results [0,0,11. For a value of ka = 1 
asymmetry of forward and backward scattering appears according to the Mie- effect. The 
ratio of forward to backward scattering Ijorw/hack increases strongly with increasing size 
parameter ka beginning at a value of ka ^ 1. "New" maxima grow out in the backward 
direction and move to the forward direction with increasing ka. The 0-dependence is not 
as striking and not as complicated as the ^-dependence since we have a superposition of the 
form A {6) -cos^ (p + B (6') -sin^ , the scattering-behavior is dominated by the strong increase 
in forward scattering described by the ^-dependence. The fact that the 0- dependence is 
described by the intervall (0, vr) is most important for the linear case and harmonics with 
odd order and differs from SHG and harmonics or even order where the 0-dependence is 
fully described by the intervall (0, |). 

The polar-plots in the case of THG are quite similar to the linear case up to values of 
ka ~ 2. The differences between Figs. ^ (b) and ^ (b) with ka = 1 reflect the stronger 
increase in the ratio of forward to backward scattering in THG with increasing ka. The 
plots in Figs. ^ (c) and ^ (c) with ka = 2 are very similar apart from one more maximum 
appearing in the third harmonic case for 6 ~ fracK2. For ka = 5 the intensities parallel 
to the direction of polarization are much larger in third harmonic than in the linear case. 
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Note, the different scales of the axes for different ka. The terms for m = 3 in THG are 
neghgble compared to the terms with m = 1. So the differences in the magnitudes of the 
hnear and third harmonic intensities are caused by the coefficients (/, 1) , afj {1, 1) and 
(1,1) only. 

The angular dependence of the intensities in second harmonic is very different from the 
linear and third harmonic cases. The disappearing direct forward and backward intensities 
and the cos^ (20) and sin^ (20) behavior produce the club-shaped structure, which is shown 
in Figure |. But the main features of the linear and third harmonic plots appear also in 
second harmonic. The plots become asymetrical with respect to 6 in the range of /ca ~ 1. 
The ratio of forward to backward scattering increases with ka and is between the linear 
and third harmonic one. Values of the forward to backward scattering ratio with different 
/ca-values are listed in Tab. II. In general harmonics with even order will show an angular 
dependence like SHG, because the Legendre polynomials of m 7^ 1 vanish at = 0, vr and 
will appear in the cosine and sine-terms in connection with n = 0,2,.,2p (p integer). 
Analogous harmonics with odd order will behave similar to the linear case. 

B. Polarization 

To learn more about the angular dependence, especially the fca-dependence of the in- 
tensity maxima, we calculate the degree of polarization using Eq. (^Ol). First, we plot the 
polarization in the case of Rayleigh-scattering (Fig. |^ (a)). For convenience and in agreement 
with the preceding section each of the figures in Fig. ^ contain the linear and third harmonic 
curves. The plots of Rayleigh- scattering (Fig. ^ (a)) are identical. With increasing size 
parameter the maximum of the polarization moves to smaller angles (the forward direction) 
in agreement with the intensity maxima. In third harmonic generation they move faster 
than in the linear case (Fig. ^ (b)). For ka = 2 (Fig. ^ (c)) there is one more maximum 
in THG in agreement with the Figs. |^ (c) and ^ (c). Increasing the size parameter up to 
5 destroys the correlation between the peaks in the polarization plots and the polar plots 
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in THG. For example the "double" peaks in the third harmonic plot with ka — 5 cannot 
be identified with "double" peaks in the polar plots but they refiect that the intensities at 
= are comparable to those at = | in contradiction to the linear case. In the linear case 
the polarization even for values higher than ka — 5 reproduces the position of the peaks in 
the polar plots and is mainly perpendicular {P (0) > 0). This is an artifact of the chosen 
refractive index and is not characteristic for linear Mie- scattering in general. For imaginary 
parts of the refractive index close to zero, the polarization is mainly perpendicular. 

The polarization plots in second harmonic generation refiect the different shape of the 
polar plots. Since the P {9 = 0) and P {6 = tt) values do not exist and the limits lim^j^o P 
and Imie^T^ P (9) are different, the plot in the Rayleigh-range is asymmetric. The decrease 
of \P\ up to zero with increasing ka is characteristic for SHG. The square of the absolute 
value of the 6- component of the electric field \Eg\ is identically zero for = 0,7r . Thus 
P (9) is a measure of the importance of the (m = 2)-term in {9, 0)|^ with respect to the 
(m — 0)-term. Thus the (m = 0)-term can be neglected with increasing ka. Only for very 
small ka, the structure in the polarization plots correlates with the shape of the polar plots. 
The first "new" maximum appears as zero in the polarization. For larger ka, however, there 
is no correlation any more. 
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C. Forward vs. backward scattering 



By computing R defined in Eq. (^) which is a measure of the difference between the 
forward and backward intensities as a function of the size parameter ka and the real or 
imaginary part of the refractive index N = n + i ■ fi' , we want to study the development of 
the asymmetry of forward to backward intensities. Of particular interest will be enhanced 
backward scattering. 

Figure ^ (a) shows the increase of forward scattering with the size of the sphere in lin- 
ear scattering. The oscillations correlate with intensity maxima in the backward direction 
resulting from maxima of the coefficients, see Probert- Jones Enhanced backward scat- 
tering appears only in the small range of 0.5 < A;a < 1 in the case of iron. The curve for 
water droplets shows no enhanced backscattering, but the overall behavior is the same as for 
metals. Varying the imaginary part of the refractive index diminishes the forward scattering. 
In the case of ka = 1 and n' > 5 an enhancement of the backward scattering occurs. The n 
dependence of R is similar. 

In the higher harmonic case the Mie- effect is strongly enhanced. The limit of one is ob- 
tained earlier than in the linear case. The higher the order, the stronger is the enhancement 
of the forward scattering. Regions of enhanced backward scattering are hard to find. In the 
case of second harmonic generation, they exist only for small n or n' around 5 and small 
ka whereas we could not find enhanced backward scattering for metal- clusters (n' ^ 0) in 
third harmonic generation so far. 
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D. Integrated intensities 



Mie-Resonances appear with increasing radius a of the sphere or absolute value of the 
wave vector k of the incident light. In this section we will only deal with the size-dependent 
resonances. 

In Fig. |I| we show the scattering efficiency Qsc of water droplets as a function of the size 
parameter. The main features are the interference structure built up by interferences between 
the incident wave and forward scattered light and the ripple structure reflecting resonant 
surface modes. Furthermore, the results suggest the optical paradoxon limfca_>oo Qext ika) = 
2. The ripple structure correlates with resonances in the coefficients (/, m) and af:^ (/, m). 



They are resonant if their imaginary part is zero. Fig. |l^ shows the first resonance of 
ag (13, 1) for a real refractive index with n = 1.5 and n' = 0. The resonance of (13, 1) 
occurs for /ca ~ 11 (see Chylek 0). For large size parameters the distance of the ripples 
can be expressed directly by the refractive index. Finite imaginary parts of the refractive 
index would damp the resonance. Then the values of the coefficients at the resonant point 



would be smaller than 1. Even if we take n' ^ 0.1 the ripple structure in Fig. |T2| does not 
appear. Accordingly, there is no ripple structure in the scattering efficency of Fe and Ni as 
a function of size as shown in Fig. |13| (a). 

In higher harmonics no structure correlated with resonances of the coefficient can be 
detected even for vanishing imaginary part of the refractive index up to a numerical accuracy 
of 10"^''. Each coefficient a^*^ is a combination of all linear electric multipole coefficients 
determined by Eqs. ( pTj) and (|12D, respectively. Correspondingly, the size dependence of 
the coefficients A^^\l,m) has many small peaks (Fig. |12| (b) and (c)). In contrast, the 
size resonances caused by multipole-combinations, and known as interference structure in 
linear scattering, are in the case of metals more pronounced in SHG (compare Fig. |1^ (a) 



and (b)). In THG (Fig. (c)) the first peak is even more dominant and the resonances 
are visible only for enhanced resolution. In both cases, the position of the resonances has a 
weak dependence on the refractive index, but the number of peaks in SHG is twice the peak- 
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number in THG ii < ka < 10. Decreasing the imaginary part of the refractive index down 
to zero changes the behavior drastically. Instead of a dominating first peak the scattering 
efficiency now shows an continuous increase for < ka < 20 with oscillations stronger than 
in the case of n' > 1. In all three cases the absolute values grow with the absolut value of 
the refractive index if metals are considered. 
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IV. SUMMARY 



By extending the classical model introduced by Ostling et al. [?] we calculated the 
angular dependence of the second and third harmonic intensities radiated by spherical metal 
clusters. Therein the n-th power of the surface charge density induced by linear polarized 
light is used as the source of the fields radiated in higher harmonics |TD|. The source 
represents the discontinuity of the electrical displacement at the surface of the sphere. 

We find that the forward Mie-scattering is in the nonlinear case even more strongly 
enhanced (see Table II). The nonlinear optical response yields a stronger size sensitivity. 
Higher multipoles contribute already at smaller size parameters ka. For example, we find 
that the light intensities perpendicular to the forward direction divided by the light intensi- 
ties in forward direction is much larger in the third harmonic than in the linear case. Since 
the values of the Legendre polynomials with azimuthal quantum numbers m = 0, ±2 vanish 
for 6* = 0, vr and due to the 0-dependence of the form A ■ cos^ (20) + B ■ sin^ (20), see Eq. 
(p!8|), the angular dependence of the second harmonic intensities is very different from the 
linear and third harmonic distributions. Especially, direct forward and backward scattering 
vanishes. In contrast, the linear Mze-scattering results (terms with m = ±1 only) and THG 
(terms with m = ±1,±3) are similar. In particular the terms with m = ±3 in THG and 
m = in SHG can be neglected, since their absolute values are much smaller than the other 
contributions. The more pronounced Mie- effect prevents enhanced backward scattering in 
higher harmonics (Figs. P). 

In order to compare with our theory an experimental investigation of the nonlinear 
Mie-scattering would be interesting. Mze-resonances play an important role in the field of 



photonic bandstructure in high refractive materials (see John ITTI), but they are studied so 
far in the linear case only. Also applications of the theory to the study of fullerenes and 
problems in biophysics (detection and growth modes of tumor cells) would be interesting. 

We will extend the theory to ellipsoidal objects to solve the classical problem of a one 
to one correspondence of scattering profile and particle shape. Furthermore, the expected 
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curvature sensitivity of the nonlinear optical response even to particles with sizes much 
smaller than the wavelength , should be detectable in the nonlinear scattering. Thus, the 
higher harmonics are a particularly useful probe for detecting small particle sizes and shapes. 
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FIGURES 

FIG. 1. The linear extinction Qext (ka) as a function of the size parameter ka and varying the 
radius of the sphere for water droplets with complex refractive index N = 1.33 + i ■ 10~^. 

FIG. 2. Difference AQ = Qext — ^- (lEef + l-E^P) / tt (ka)'^ of the extinction efficiency and 
the modified forward scattering amplitude, as a function of the size parameter ka. This illustrates 
how our numerical calculations fulfill the optical theorem. 

FIG. 3. Polar plots of the linear scattered intensities for Fe: (a) scattered linear intensity 

{9, (p) in units of IQ-^^ with ka = 0.001, (b) {6, (p) in units of lO"'^ with ka = 1, (c) {9, (p) in 
units of 1 with ka = 2, and (d) {9, <j)) with fea = 5 in units of 1. The inset shows the scattering 
geometry. The direction of the incident light is determined by = with polarization along the 
axis defined by = and 9 = 7r/2. 

FIG. 4. Polar plots of the third harmonic (THG) intensities for Fe: (a) scattered 3a;-intensity 
huj {0, 0) in umits of 10~^^ with ka = 0.001, (b) Isw {9, (p) in units of 1 with ka = 1, (c) Is^^ {9, cf)) 
in units of 10~^ with ka = 2, and (d) I,^ {9, 4>) in units of 10~^ with ka = 5. 

FIG. 5. Polar plots of the second harmonic (SHG) intensities for Fe: (a) scattered 2a;-intensity 
hu {9, 4)) in units of 10"^^ with ka = 0.001, (b) huj {9, 4>) in units of 1 with ka = 1, (c) (9, 4>) 
in units of 10^ with ka = 2 and, (d) I^j {9, (p) in units of 10^ with ka = 5. 

FIG. 6. Degree of polarization P for the linear and third harmonic case as a function of the 
polar angle 9 (a) for ka = 0.001, (b) for ka = 1., (c) for ka = 2, and (d) for ka = 5. 

FIG. 7. Degree of polarization P for the second harmonic case as a function of the polar angle 
9 (a) for ka = 0.001, (b) for ka = 1, (c) for ka = 2, and (d) for ka = 5. 

FIG. 8. Degree of Mie- asymmetry R (defined as R = {Iforw ~ hack) / {^forw + hack)^ see 
Eq.(23)) in the linear case (a) as a function of the size parameter ka for refractive indices re- 
ferring to Fe, Ni and water droplets, (b) as a function of the real part of the refractive index n, 
and (c) as a function of the imaginary part of the refractive index n' for size paramters ka = 1 and 
ka = 5 using values of fl, n' for Fe and Ni. 
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FIG. 9. Degree of Mie- asymmetry R as a function of the size parameter ka for (a) SHG, and 
(b) THG with the refractive indices referring to Fe and Ni, respectively. 

FIG. 10. Degree of Mie-asymmetry R as a function of the real part of the refractive index h. 
For ka = 1 and ka = 5 for (a) SHG, and (b) THG. n' corresponds to Fe and Ni. 

FIG. 11. Degree of Mie-asymmetry Ras a function of the imaginary part of the refractive index 
n'. For ka = 1 and ka = 5 for (a) SHG, and (b) THG. n corresponds to Fe and Ni. 

FIG. 12. (a) Real and imaginary part of the expansion coefficient of the scattered linear fields 
a^' (/ = 13, rn =1) as a function of size parameter ka and for a real refractive index h = 1.5. 
The first peak represents the lowest resonance of the multipole mode with I = 13. The real and 
imaginary parts of the expansion coefficients of the SHG and THG fields and , respectively, 
as a function of the size parameter ka are shown in Figs, (b) and (c). The results in (b) and (c) show 
that the resonance peaks of the linear coefficients contribute to the higher harmonic coefficients: 
The maxima and minima of A^^^ and A^^^ are essentially determined by the resonances of the linear 
coefficients a^. 

FIG. 13. Scattering efficiencies Qsc as a function of size parameter ka for (a) linear scattering 
(b) SHG, and (c) THG in the case of Fe and Ni, respectively. 
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TABLE 1. 
and THG. 



TABLES 

Ratio of forward to backward intensities as a function of ka in the linear case, SHG 





I forw / 1 


back 


ka 


lin 


SHG 


THG 


1 


1.4 


2.7 


15 


2.1 


8.7 


10 


113 


4.8 


64 


136 


206 



TABLE IL Complex refractive indices N = n + in' for iron and nickel used in this work at 
u, 2u and 3a;. The wavelength is A = 617 nm and a; = 27r/A • c. 





iron 


nickel 




n 


h' 


n 


n' 




2.88 


3.05 


1.99 


4.02 




1.69 


2.06 


2.01 


2.18 


?>UJ 


1.49 


1.41 


1.29 


1.89 
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